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Abstract 



Diffusion on a T fractal lattice under the influence of topological biasing fields is studied by 
finite size scaling methods. This allows to avoid proliferation and singularities which would 
arise in a renormalization group approach on infinite system as a consequence of logarithmic 
diffusion. Within the scheme, logarithmic diffusion is proved on the basis of an analysis of various 
temporal scales such as first passage time moments and survival probability characteristic time. 
This confirms and puts on firmer basis previous renormalization group results. A careful study 
of the asymptotic occupation probabilities of different kinds of lattice points allows to elucidate 
the mechanism of trapping into dangling ends, which is responsible of the logarithmic time 
dependence of average displacement. 
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1 Introduction 

A particle diffusing under the influence of a random biasing field can exhibit slow diffusion, with average 
displacement, R, growing in time, t, as some power of Int [1-7]. A remarkable exact demonstration of 
logarithmic diffusion for a particle hopping on a one-dimensional chain and subject at each site to an 
independent random bias is due to Sinai [7]. There is also evidence, based upon numerical results, that 
logarithmic diffusion replaces the anomalous power law one, when particles diffuse on fractal structures, 
like the incipient infinite cluster of percolation, and are subject to a uniform biasing field [8,9]. In this 
case, though the field is not random, the complexity of the fractal structure conspires with the former in 
determining possible localization effects, e.g., by pushing and trapping the moving particle into dangling 
ends. So far, this phenomenon has been modelized, in most cases, by diffusion along comblike structures 
with teeth of variable length, where a bias pushes the particle towards the tips [10]. Despite the efforts 
produced, a comprehensive and satisfactory mathematical description and a general physical understand- 
ing of the phenomenon of logarithmic diffusion arc still lacking. Present insight is mostly limited to the 
example of a particle diffusing in ID and subject to a random bias [11]. The case of diffusion on fractal 
structures under constant bias, which is our concern here, is more problematic. Comblike models of 
this latter type of diffusion are again basically one- dimensional, and rather far from the complexity of 
self-similar structures. The accepted picture of bias induced logarithmic diffusion on a fractal, through 
the mechanism of trapping into dangling ends, is intuitively appealing but rather vague. So, it is worth 
to further elucidate and test it, in order to have a more reliable description of such a process. 
Recently a study of such trapping mechanism has been accomplished by a dynamical renormalization 
group (RG) analysis of diffusion on deterministic fractal lattices under the influence of a topological bias 
[12]. Deterministic self-similar structures, though being a mathematical schematization, retain various 
fundamental features of real random systems, with the advantage of allowing analytical solutions of ran- 
dom walk processes [13-15]. Furthermore, the RG has demonstrated to be a powerful tool in treating 
dynamics on self-similar structures [14,15], also in the very peculiar singular conditions implied by log- 
arithmic diffusion [12]. Indeed, one should notice that logarithmic time behaviour implies a singular 
structure for the RG transformation. A standard dynamical RG yields the time rescaling factor P, under 
length rescaling I, which is consistent with power law anomalous diffusion, R ^ i^l^' . So, logarithmic 
diffusion clearly corresponds to the z +oo limit, causing ^ oo as well, and thus a singularity in the 
RG transformation. As claimed in ref.[12], the control of a singular RG transformation should be the key 
to extract the correct logarithmic time dependence of R. On the other end, a singular RG transformation 
can be useful only if the RG is carried out exactly [12]. Such a requirement and the consequent need 
of handling extremely strong singularities reveals a very hard task, implying, by necessity, some ad hoc 
simpliflcations and assumptions in the RG approach [12]. 
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In alternative to RG, finite size scaling (FSS) analysis of non-equilibrium dynamical processes [16] can be 
exploited for the study of power-law anomalous diffusion. Recently a remarkably efficient and straight- 
forward approach to diffusion on finite self-similar trees and hierarchical combs has been proposed [17]. 
In this approach, through an exact decimation of the master equations and consequent rescalings, the 
FSS of first-passage time (FPT) moments and survival probability (SP) can be studied. 
A FSS based phenomenological RG study of anomalous diffusion offers several advantages with respect 
to conventional RG strategies. In particular one does not need to care about proliferation problems in 
the dynamical equations. Thus, one can legitimately wonder whether a similar approach could be used 
to confirm and elucidate logarithmic anomalous diffusion on fractal structures of the type considered in 
ref.[12]. 

In the present paper we apply the methods of ref.[17] to address diffusion on deterministic T-like finite 
fractal trees in two dimensions, when a topological bias field is acting. We thus generalize these finite 
size scaling methods [17] to a system exhibiting singular dynamical scaling. At a first sight, this singular 
scaling could seem to prevent any effective analysis altogether. To the contrary, we are able to derive 
the FSS behaviour of FPT moments, SP, first-passage probability (FPP) and probability of being at a 
certain site at time t (the site occupation probability). All these quantities turn out to be controlled by 
the same time scale, t* ^ e"^, L being the linear size of the finite lattice. This confirms that, in the 
infinite system, diffusion is logarithmically localized within a distance R'^lnt, after a time t, as already 
conjectured [12]. Furthermore, as discussed below, the existence of a single time scale, t*, characterizing 
both the fixed observation-point ensemble (yielding the average over the first passage time at a fixed 
given site, i.e. the FPT moments) and the fixed observation-time one (concerning the average over the 
walker displacements at time t, {R{t)) ), supplies concrete support to the intuitive hypothesis that bias 
induced trapping into dangling ends is directly responsible of logarithmically slow diffusion on self-similar 
structures. 

As a further point of interest, our results show that the presence of a singularity in RG transformations 
at the fixed point does not prevent the FSS hypothesis to hold in this case. This was not a priori obvious 
in such a circumstance. 

As anticipated above, from a methodological viewpoint, the strategy of studying logarithmic diffusion by 
FSS allows to avoid some of the difhculties involved in renormalising the infinite model. Since finite size 
analysis requires to consider the infinite time limit of physical quantities while keeping linear dimensions 
finite, asymptotic time behaviours are extracted while divergences, induced by the singular time rescaling, 
still have to develop fully. Hence no ad hoc assumptions are needed to handle strong divergences, as it is 
the case in infinite systems [12]. 

This article is organised as follows. In the next section we introduce the model of biased diffusion on 
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deterministic T fractal, and discuss quantities relevant to the description of the random walk process. In 
the third section the decimation technique of the finite set of master equations is developed and applied 
to our model. Phenomenological RG transformations of FPP and SP are then derived and the scalings 
of FPT and characteristic time are calculated. The methods of section 3 are reconsidered in section 4 
and used to extract the long time behaviour of occupation probabilities. Several variants of the model, 
differing in initial and/or boundary conditions arc also examined in section 4, to show that these differ- 
ent conditions do not influence the main dynamical scaling properties. Finally section 5 is devoted to 
discussion and conclusions. 

2 Deterministic fractal model of biased diffusion 

Diffusion on random ramified self-similar structures can be fruitfully investigated by the study of random 
walks on deterministic fractal lattices [13-15]. Often these offer the advantage of being simple enough 
to allow analytical treatment. For definiteness we will focus on a particular realization of our model, 

the so-called T fractal (see Fig. 2.1), a hierarchical tree in 2D, whose fractal dimension is J = In 3/ In 2 
[18]. On such a lattice at some finite generation order n, a particle diffuses by hopping between nearest 
neighbour (nn) sites. The probability, Pi{t), for the particle to be at site i at time t, obeys the master 
equation: 

m + r) = Pi{t) + J2 [WijPiit) - Wj,Pi{t)] , (2.1) 

where the sum on the right hand side nms over nn, j. of site «, and Wij is the probability that a 
particle hops from j to i in time t. Here it is assumed that Wij = W+ or W-, according to whether 
one goes from site j to i following the bond arrow or not, respectively. Arrows are all pointing in the 
direction of increasing "chemical" distance from the lattice's leftmost site (site 1 of Fig.2.1). Thus the 
following consistency relations have to be satisfied: 2W+ + W- < 1, W± < 1. Moreover, we assume 
W = W-/W+ < 1, so that the bias tends to drive the particle towards longer chemical distance from 
site 1. This is an example of topological bias field, namely a field having a particular direction in the 
topological space where the hopping process takes place. This is different from the so called Euclidean 
bias, which is uniformly directed in the Euclidean space onto which the lattice in embedded [13]. 
We choose site 1 to be the starting point for particle motion: so, the initial condition of eq.(2.1) is 
^4(0) = 6i_i. We also impose the rightmost end-tip site of the lattice, is, (e.g. number 4 or number 10 of 
Fig. 2. la or b, respectively) to be a fully absorbing sink, i.e. a position from which the particle can not 
hop back to other sites of the T fractal {Wji^ = Vj). 

Due to the absorbing site, the SP has an exponential decay whose characteristic time is related to the 
linear size of the lattice. At generation n, a T fractal has linear size L = 2"+^. In standard anomalous 
diffusion processes, where the average displacement, R, scales as ~ t^^^, the characteristic time, t*, is 
expected to scale with L ast* L^. Relying on previous results for the same model without sink and in 
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the infinite size hmit [12], we expect in the present case a logarithmic locaHzation of motion too, thus a 
scahng t* e"^, with a > 0. 

Moreover, since is is absorbing, its occupation probabiHty, Pi^ (t) , coincides exactly with the first-passage 
probability, Pppit), at time t for that site. Pi^ is also closely related to the survival probability, S{t), 
namely the probability that the particle has not yet reached the output (sink) site at time t. Indeed, 



S{t) = 1 - XI PFp{kT) , (with t = mr). (2.2) 

k=0 

Putting T = 1, and applying the time discrete Laplace transform [14,19], 

oo 

P°(a;o)= XP.(m)(l+c.o)-^-" , (2.3) 

m=0 

to eq.(2.1), we get: 



^*H = E^^:'-M + ^'.i ' (2-4) 



where oj = coo/W+, P^{coo) = Pk{co)/W+, and the Kronecker's delta reflects the initial condition Pi{0) = 
6i^i. According to eq.(2.1) one should have a{i,co) = w Mi, in eq.(2.4). However, we already introduce 
these coefficients in view of the fact that sites with different coordination are inequivalent under lattice 
decimation and consequent rescahng, as we shall see later. It will turn out that a{i,u)) = Qi(ijj), ay,{ijj) 
or ao(a;), according to whether i has coordination 1, 3 or coincides with site 1, respectively. The only 
parameters entering eq.(2.4) will thus be a = (aoj 0^1,0:3)5 and W. 

When the sink site is is considered, we see immediately that Laplace-transforming Pi^{m) = Pppim), 
provides the generating function for the FPP, indeed: 

{coo + l)^-jP^\ =- V(m + l)Pfp(m) = -(m)-l , (2.5) 

(m) = {t) is the average FPT. Higher moments of FPT come from higher derivatives of Ppp- E.g., by 
deriving Ppp two times one finds: 

(wo + 1)^ (wo + l)^P^p(wo) I = V (m + l)(m + l)Ppp{m) = (m^) + 2(m) + 1 . (2.6) 

^ ' m=0 

Clearly most informations about the diffusion process are contained in PP(a;o), or equivalently in Pi(u)). 
Thus, the analysis of these quantities will be our main concern below. 

3 First-passage time and characteristic time scaling 

The first-passage probability for the T fractal with an absorbing point at the right end tip, can be 
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calculated by successive decimations of the Laplace-transformed master equation. Consider the zeroth- 
order tree of Fig. 2. la, with site 1 being the input site and site 4 the output one. The explicit form of 
eq.(2.4) reads: 



Solving this system for P4 yields: 



f {ao{oj) + l)Pi = WP2 + 1 

(as{uj) + 2 + W)p2 = Pi + Wp3 
{ai{uj) + W)p3 = P2 

[UJ+^)P, = P2 

ai{Lj) + W 



(3.1) 



(w + ^) {(ao(w) + 1) [(a3(w) +2 + W) (ai(w) + W)-W]- (ai(a;) + W) W} 
Furthermore, the Laplace-transformed master equation system for the first order tree of Fig. 2. lb is: 

' (aoiuj) + l)Pi = Wp2 + 1 
{ai{u) + W)Pi=Pj 



(3.2) 



ia^ito) + 2 + W)Pi = Pk + Wj2Pj 

(aaiuj) + 2 + W)p8 = P4 + WPg 



(3.3) 



{uJ + ^)P,o=Ps 



Decimating the last generated sites, namely 2, 3, 5, 6, 8, and 9, and eliminating the respective p's from 
system (3.3), we go back to a lattice which is equivalent to the first order one, after rescaling lengths by 
a factor 2. After suitable redefinitions and site index reordering, the system for the surviving P's can be 
reduced to a form analogous to (3.1), 

f (4')(a;) + l)Pi(') = VK(i)pi'' + l 

(41) (^) + 2 + 1^(1) )p(i) = p(i) + W^'^P^''> 

(a«(a;) + W^W)Pi^)=pf^ 

with the following identifications: 

Ai)( ,^ (ao(u^) + 1) [(«i(u^) + W){a:,{u;) + 2 + W)-W] 
a}, {uj) = — 1 — W 



(3.4) 



a^lHuj) = [aiiu) + W){aj.{uj) + 2 + W)-2W -W^ 

^(1) . ^ ^ (M^) + 2 + W) [{a.ju) + W){as{io) +2 + W)-W] _^_^^_^2 
^ ^ ' ai{u) + W 



(3.5a) 
(3.56) 
(3.5c) 
(3.5rf) 



Q!l(w) + W 



(3.5e) 



iai{uj) + W)ia3{uj) + 2 + W)-W 
with i' = 1 if i = 1, z' = 2 if i = 4, and so on. The apex (1) in front of all left hand quantities indicates 
that they have been rescaled once after a decimation. A few remarks about these RG transformations: 
in first place, quantities relative to the sink site do not rescale, as it appears from the last of eqs.(3.4). 
In addition one should notice that different a-^^'s undergo different transformations, the dependence on 
the original aj's being nonlinear. This is a standard memory effect associated with dynamical coarse 
graining. Now the necessity of introducing different a-coefHcients becomes clear. 

Starting with an nth-order lattice, we perform n successive decimations and rescalings in such a way to 
arrive to the zeroth-order system. Then, solving the latter for P4, the FPP generating function, we find: 



(3.6) 



Moments of the first passage time arc obtained by the scries expansion of Pfp- Indeed assuming that 
Ppp is expandable around co = (which is justified for a finite system) i.e. 



we see immediately from (2.5-6) that pi = — (t) — 1 and p2 = {t^)+2{t) + 1 (notice that, since 

1 dpFpiuJo/W+) dpFp{uj) 



(3.7) 

dujo 



^ ) Pi ^^c related to FPT moments without any further rescaling). Thus, 
once the relation between the pi's of two successive T Fractal orders is known, one can derive the scaling 

of (t). Since we arc interested only in the average FPT we focus on pi rescaling. To this aim we have to 
expand the right hand side of (3.6) and inherent quantities up to first order in w about lu — 0. 

To calculate a("^(tj) wc first consider its series expansion in powers of a; about uj = 0. Writing a["\uj) = 

00 

^^a^"^''^a;'^, and considering (3.5b-d), we are able to derive the relation between and a("~i)'i, 

fe 

which are needed to calculate the scaling of pi. By keeping everywhere only leading order in 
(W^n) ^ 1^2" ^ w^/^), in the explicit relations, we get: 

(1 + M^)4"-')'' _ (l + iy)2"-i 



(n),l ^ 



a 



(n),l 



^(n),l ^ 



(1 + W)2aJ""^^'^ ^ (1 + H^)2" 
(1 + W)24""^^'^ _ {1 + W)2" 



(3.8) 



Disregarded terms in (3.8) are o{W ), which is negligible for n » 1, since W < 1 in the biased case. 
Expanding (3.6) consistently with the approximation leading to (3.8) then yields: 

W+ 



I^(") - wa^"^'' + 0{uj^) 
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(3.9) 



So {t) scales as 

(i)^"^ + l = ^j^((i)^"-^^ + l) , (3.10) 

and finally we have: 

The survival probability, S{m), defined in (2.2) decays exponentially in m as S{m) oc exp(— m/f*), 
at long times, for any finite nth-order T fractal, due to the absorbing boundary conditions. So, we 

want to investigate the relation between the characteristic time, t*, and the FPT. Ppp{m) decays also 
exponentially, because it is just the discrete derivative of S(m): Pppim)'" cx " (e^/'" - l) e"™/*". This 
implies that Pfp{u)) has a simple pole in ujc = -tt— (e"^/*" — 1) in the limit w — » 0, which corresponds 
to the long time behaviour of Pppim). 

The characteristic time grows with the T fractal size n, f*"^^^ + 00 ; hence u:c—^^~ and can be located 
by looking for the simple pole of Ppp in the limit a; ^ [20]. Expanding numerator and denominator 
up to order w on the r.h.s of (3.6), and exploiting (3.8) for the expansion of a^"^(a;), we eventually arrive 
to: 



Ppp{^)^- , ^ . (3.12) 



The first simple pole = is not related to long time behaviour, since — < — 1. The second 

vf(") T4^2" »-+oo -1 

simple pole is Wc — r^rr — ; 777^ — = —-, ttttt- Since — ttt , then 

^ ^ 2oif>^^ 2(1 + VF)2" {\ + W)L W+t*' 

2"+i Ij 

So diffusion has a characteristic time scale growing exponentially with the system size, and the dynamical 
exponent, z, diverges as expected. Moreover scales like and, consistently with infinite lattice 

result for the average displacement [12], the typical distance, L, travelled in time t is i ~ Int. 
All the above results provide evidence that dynamics is controlled by a single time scale determined by 
the trapping mechanism into dangling ends. In truly infinite fractal structures [12] dangling ends of all 
length scales exist, thus dynamics is logarithmically localized. 

The existence of a single characteristic time suggests also that fixed time and fixed point of observation 
ensembles are substantially equivalent in this model, in accordance to what happens in most models of 
power-law anomalous diffusion on self-similar lattices [13-17]. 

Finally, two important remarks concerning the approximate calculations are in order. In first place, due 
to approximation to o(W^^" ^) introduced to obtain eqs.(3.8) and consequently (3.9) and (3.12), results 
concerning (f) and t* can not be extended to the region < 1, of cross-over to the standard anomalous 
diffusion regime corresponding \,oW = \ [12,14]. Secondly, we assumed Pfp{co) to be expandable about 



a; = to get (3.9), while we have then seen it has a pole Wc"— — ?°0~. However, since Wc is finite for 
any finite n, no matter how large, this will cause no inconsistency, because for oj close enough to 0, an 
expansion will always be allowed. On the other end, this feature implies that the limits t — > +oo and 
n +00 (i.e. L oo) can not commute, because when L ^ oo a singularity of Ppp sets exactly at 
Lu = 0, thus results (3.9-11) do not extend to that limit. This problem does not exist in principle for 
the asymptotic analysis of FPP: indeed, the derivation of (3.12,13) involves no analyticity assumption 
for Ppp- However, the determination of the n — > +oo limit remains a serious task in its own due to the 
presence of strongly divergent terms in RG equations [12]. 

4 Long time behaviour of the occupation probabilities and trapping mechanism 

The knowledge of the occupation probability, Pj{t), of a generic site j of the T fractal at its nth-order, at 
time t, supplies further information about the dynamics. Pj (t) is also expected to decay exponentially at 
long times, and we are interested in its characteristic time, which should scale like t* . Furthermore, since 
sites with different coordination are dynamically inequivalent (look at different a^'s rescalings in (3.5b- 
d)), such a feature should extend, in some way, to their occupation probabilities. When these show an 
asymptotic exponential decay, they should also have an amplitude factor dependent on site coordination 
and position, and on the bias parameter W as well. The nature of such dependence may provide further 
insight into the diffusion mechanism. 

Asymptotic analysis applied to Pppit) can be extended, with moderate additional effort, to get also the 
asymptotics of a generic Pi{t). 

Consider again the nth-order T fractal with initial condition at site 1 and absorbing boundary condition at 
the rightmost end-tip site. We perform n successive decimations of the master equations (2.4), eventually 
arriving to a system, like (3.4), for the surviving p(")'s, rescaled n— times according to (3.5). Now the 
system can be completely solved for the occupation probabilities of surviving sites, yielding: 



^3"'H = ^ (4.1a) 
( a^"^ (w) + ( af^ (w) + 2 + W^^A - 



1 



J-i"^(^)= ^ ((^ + H/ ^^^^^^ 
A(") = (a^"^(a;) + 1) [(a^"^(w) + W("))(4"^(a;) + 2 + W^")) - W^")] - W(")(aJ"^(w) + W^(")) 
and ^4(0;) as in eq.(3.6). Using (3.5e), we then recover the original, unrescaled, P's as: 

" [a''^^ {uj) + W^^^) (af ) (w) + 2 + W^^^ ) - W^^^ ' 



PM = n 



fe=0 



a^'=^(a;) + VF(fe) 
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where i' is the index that site i takes after n decimations. Again P{u}) is a fraction and, being concerned 
with long time behaviours, we expand both numerator and denominator to leading order in uj. In this way 
we expect to find the single pole corresponding to the characteristic time, like already seen for Ppp{u)). 
As additional task here, the product preceding P^") on the r.h.s. of (4.2), has to be expanded. To leading 
order in w, this can be accomplished using the series expansions of a'^'^^{uj) about w = and exploiting 
equations (3.8) for the coefficients a\'^^'^- Relations (3.8) are drawn within an accuracy of o{W^"~ ), 
which thus applies also to the results below. We shall omit the details of such cumbersome calculations, 
and give directly the asymptotic expression of occupation probabilities: 



(1 + W) 2"+i 



For purposes which will be clear later on, sites have been labelled with the indexes they have after n — 1 
decimations, consistently with Fig. 2. lb. So 1 is the injection site, 10 is the absorbing one, 7 is at the tip of 
the central branch and 4 is at the lattice center (they correspond to sites 1,4,3,2 of Fig. 2. la, respectively). 
Note that the characteristic time, , is the same we found in (3.13), for every site, while the amplitude 

dependence on W is rather varying from site to site. 

In order to have further terms of comparison it is worth deriving the asymptotic P{t) form for a few 
other sites too, namely those which disappear due to the nth decimation, but have different chemical 
distances from both origin and sink. For definiteness we focus on sites 3,6 and 9 which are tips, and on 8 
which has three nn (indexes refer always to Fig. 2. lb). We suppose to arrest decimation at the (n — l)th 
stage, obtaining a lattice configuration like Fig. 2. lb. The site occupation probabilities to be calculated 
can be drawn by solving the corresponding (n — 1)— times decimated master equation set for Pj;"~^\u;), 
and using also results (4.1a-c). For example, site 6 is at the same chemical distance from 1 as 7, then 
6 and 7 are completely equivalent and their P(t)'s show identical long t behaviour. On the other end , 
site 9, though equally far from the origin as 7, is influenced by the absorbing boundary condition on the 
branch it belongs to and exhibits different features. Indeed: 

, . p(»-i) 
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(see e.g. system (3.3)), and since 

{a^"-^\uj) + W("-i))(4"-^^(a;) + 2 + VF("-i)) - 



n = .,„- ir^il"'^'"" pt" = pt> 



we have 



n-2 



^9(0;) = n 



fe=0 



(a^'=) (w) + ) (a^''^ (w) + 2 + VF(^) ) - W^^) 



af)(a;) + W^W 



a^")(a;) + Ty(") 



(a^"-')(a;) + W^("-i)) AO 



(4.4) 



Applying tlie above outlined asymptotic analysis to Pg, one finally arrives to 

notice the W^"^ ^ extra factor, comparing to eq.(4.3c). In a similar way we calculate occupation proba- 
bilities of another tip site, number 3, which is closer to the bias origin, and another site with coordination 
3, number 8, 

^sW-H^W^f^e-A^ , (4.6a) 

The above results, (4.3,5,6), again point out how in a branch of linear size L, all occupation probabilities, 
as well as the FPP, decay with the same characteristic time t* ^ g|ii^l-'^/2^ jj^ infinite structure 
with branches of any length scale, the larger is the branch, the exponentially longer will be the time the 
particle spends in it, due to trapping effects caused by the bias. This will cause a logarithmic slowing 
down of diffusion. 

Strong evidence that the bias drives the particle into the deepest dangling ends can also be drawn by 
looking at the amplitude dependence of occupation probabilities. Indeed we notice that the P{t) of tip 
sites which are at maximum chemical distance from the bias source, have an amplitude factor proportional 
to 1/2" (see eq.(4.3c); an identical result holds for site 6). The number of such sites grows just as 2" 
with the order n of the T fractal, because the branch containing the absorbing site has not to be taken 
into accoimt due to its depleting nature (compare for example eq.(4.3c) with (4.3d) and (4.5)). Hence we 
infer that, after a transient time, probability accumulates into these deep dangling ends, redistributing 
then slowly against the biasing force to sites with higher coordination and/or in less favourable positions. 
Amplitude factors for the probability of every other kind of site are in fact proportional to W''/2", where 
the exponent r] depends manifestly, though non trivially, on the chemical distance between the site under 
consideration and those deepest dangling end, namely on the number of unfavourably biased bonds one 
has to cross going from the latter to the former. In an infinite size limit, occupation probabilities of sites 
with coordination > 1 will then be vanishingly small compared with those relative to the dangling ends, 
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in agreement with our hypothesis. 

Finally we would like to point out the peculiar features of the sink site. This site can not be considered, like 
other tips, as a probability accumulating point, of course, and the same applies to its neighbouring sites 
(e.g. number 9 in Fig. 2. lb). Such circumstance reflects on the W^" and W"^" ^ prefactors in eqs.(4.3d) 
and (4.5), respectively. Moreover Pio (i.e. Ppp) has actually two simple poles: ujd = —^^^ anda'c2 — — 
1/W+t* , as it appears from eqs.(3.6,12). Therefore two characteristic times exist, the first deriving from 
Wci, very short and connected to the transient time when the initial flux of probability entering the 
side branch is immediately absorbed by the sink. The second, , is quite long and corresponds to the 
asymptotic regime when probability formerly accumulated into deep dangling ends slowly reflows in the 
whole lattice and then is absorbed by the sink. 

To probe a possible dependence of diffusion process on initial conditions, we have studied the same model 
choosing different input sites, namely no longer site 1, but either 4 or 7 of Fig. 2. lb, successively. We 
arrived at exactly the same results for {t)^^\ f*, Pppit) and occupation probabilities. Hence we claim 
that the above results and considerations hold for our model independently of initial conditions. 
The description of the dynamic model we are concerned with, should be independent on the particular 
realization of the bias pattern. To get convinced of this, a less asymmetric situation than that presented 
in Fig.2.1, has also been considered. Again, starting from a T fractal, we suppose to orient the arrows of 
topological bias as in Fig.4.1. This new choice, together with the conventions on Wij specified in section 
2 [21], implies that the bias pushes a diffusing particle away from the T fractal central site (number 2 in 
Fig. 4. la, or 4 in Fig. 4. lb). In addition the three tip-sites (e.g. 1,7,10 in the 2-order lattice of Fig. 4. lb) are 
assumed to be absorbing, and the central one to be the input point of the walker. Master equations and 
their Laplace transforms are equivalent to (2.1) and (2.4) respectively, with the difference that now ao{co) 
refers to the central site, which has three nn, but is not dynamically equivalent to other 3-coordinated 
sites. The method outlined in section 3 and the present one applies of course to the latter model too. 
Considering an nth-order tree one finds, for relevant quantities: 

(1 + W) 2" 

^ ^ 1 - 1^2 3 2« ^ ' 

P3,6.9(0 ^ ^ ^e-*/< (4.7e) 
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Site indices, shown in Fig. 4. lb, are those which would apply after n—1 decimations. Taking into account 
that the linear size of the nth-order lattice is 2"+^, but that each of its three branches is indeed equivalent 
to the (n—1)— order lattice of Fig. 2.1, and its length is thus L = 2", then results (4.7) are totally consistent 
with previous ones. This is not surprising in view of the threefold symmetry implied by the conditions 
we imposed. 

So, our analysis applies to different realizations of the present model. It can also be shown that our 
approach can be extended to a wider family of ramified fractals, like for example that in Fig. 2. 2, in d > 2 
dimensions. 

5 Conclusions 

We have considered the problem of random walks on deterministic self-similar trees in the presence of a 
topological biasing field. Relying on FSS ideas, wo have developed a method to extract the FSS behaviour 
of the average FPT, (i), and characteristic time, t* . The method is phenomenological and is based on 
successive decimations and rescalings of the set of master equations for a finite nth-order lattice, {t) and 
t* are found to scale in the same way with the system size, L, i.e. (t) ~ ~ ' "^^u^' dynamics is 

controlled by a unique time scale growing exponentially with the tree's size. Since infinite fractal trees 
contain branches of any length scale, where the particle can be trapped in, we argue that the resulting 
slowing down should indeed lead to the logarithmic diffusion conjectured for infinite systems. Moreover, 
the existence of a unique time scale reveals a substantial equivalence between the fixed time and the fixed 
point of observation ensembles, governing the behaviour of {R{t)) and (t), respectively. The scaling of t* 
implies also a divergence of the dynamical exponent z and thus a singularity in standard RG transfor- 
mations for infinite lattice. The possibility to carry on a phenomenological RG strategy in the presence 
of singular transformations, via the expedient of considering systems of finite size, should be seen as the 
main achievement of our technique, in view of the serious complications involved in the RG analysis of 
infinite systems [12] . It is also worth noting that we produced here an interesting example of the validity 
of the FSS hypothesis, notwithstanding the existence of a singularity in the RG group. 
Finite size analysis has been pushed further here to calculate the long time behaviours of FPP and occu- 
pation probabilities, by combining asymptotic expansion of the respective Laplace transforms with RG 
rescalings. All probabilities are found to decay exponentially at long t , again with the same character- 
istic time t*, as expected. This fact, of course, supplies stronger confidence to the picture of logarithmic 
diffusion generated by the mechanism of bias induced trapping into dangling ends. Amplitude factors of 
occupation probabilities offer further insight into this dynamic process. Indeed, the occupation proba- 
bilities of sites located in the deepest dangling ends have amplitude proportional to 1/2", which is just 
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the inverse of their number in an nth order lattice. Moreover, probabiHties relative to sites with higher 
coordination or in a less favourable position with respect to the biasing field, have amplitudes propor- 
tional to /^^ iyV < 1) , with r] depending on the chemical distance between the deepest dangling 
ends and the site under consideration. Thus probability first accumulates into the longest branches of the 
fractal under the bias drift; then it very slowly reflows against the bias to the whole lattice, generating 
the logarithmic slowing down of diffusion. 

The above results have been demonstrated to hold in general for the system \mdcr discussion, no matter 
which initial and boundary conditions are chosen, and independently of the particular realization of the 
model. Our methods can be applied to a wide variety of finitely ramified fractals [18], for which they 
provide the same picture of logarithmic biased diffusion outlined here. 
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FIGURES 



FIGURE 2.1 The first three successive iterations of the T fractal, (a), (b) and (c) respectively. The arrows on 
each bond indicate the direction of the topological bias. 



FIGURE 2.2 Example of T fractal with coordination 4 or 1, in d = 2. Analogous lattices can be drawn in d > 2. 



FIGURE 4.1 The first two iterations of the symmetric T fractal lattice, (a) and (b) respectively, with a topological 
bias on each bond. 
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